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Quantum discord, a measure of genuinely quantum correlations, is generalized to continuous variable sys- 
tems. For all two-mode Gaussian states, we calculate analytically the quantum discord and a related measure of 
classical correlations, solving an optimization over all Gaussian measurements. Almost all two-mode Gaussian 
states are shown to have quantum correlations, while for separable states, the discord is smaller than unity. For 
a given amount of entanglement, it admits tight upper and lower bounds. Via a duality between entanglement 
and classical correlations, we derive a closed formula for the Gaussian entanglement of formation of all mixed 
three-mode Gaussian states whose normal mode decomposition includes two vacua. 
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Entanglement, non-classicality and non-locality are among 
the pivotal features of the quantum world. While for pure 
quantum states these concepts are all equivalent, being like 
three facets of the same gemstone, they coiTespond to differ- 
ent resources in the general case of mixed states. Namely, 
while entanglement plays a central role in bipartite and mul- 
tipartite quantum communication [IJ, its necessity for mixed- 
state quantum computation is still unclear Conversely, 
several recent studies have shown that separable (that is, not 
entangled) states, traditionally referred to as "classically cor- 
related", might retain some signatures of quantumness with 
potential operational applications for quantum technology 1 3 - 
IS). One such signature is the quantum discord |7|, which 
strives at capturing all the quantum correlations in a bipartite 
state, including - but not restricted to - entanglement. Pro- 
gresses in the demonstration of quantum protocols for the ma- 
nipulation and transmission of information, as well as in the 
quantification of entanglement, which is closely related to the 
efficiency of such protocols, have been amply recorded in the 
case of finite-dimensional systems, mainly qubits, as well as 
infinite-dimensional systems, where correlations between de- 
grees of freedom with continuous spectra are exploited [8J. 
However, there persists a fundamental gap between finite and 
infinite dimensional systems concerning the investigation of 
more general measures of quantumness versus classicality |9l- 
For Gaussian states, the workhorses of continuous variable 
quantum information, such an investigation would be espe- 
cially valuable, since in view of the positivity of their Wigner 
distribution these states have been sometimes tagged as essen- 
tially classical. 

In this Letter we endeavor to bridge this gap. We define the 
quantum discord for Gaussian states and explicitly solve the 
optimization problem involved in its definition, constrained 
to measurements that preserve the Gaussian character of the 
states. We derive a closed formula for the ensuing Gaussian 
quantum discord and for a related measure of classical corre- 
lations fTOl on all two-mode Gaussian states. We prove that 
these quantum correlations are limited for separable Gaussian 
states, yet they are nonzero for all but product states. For en- 
tangled states, quantum discord is found to admit tight upper 



and lower bounds, given by functions of the Gaussian entan- 
glement of formation [11]. Exploiting a duality between en- 
tanglement and classical coiTelations |12|, we derive an an- 
alytical formula for the Gaussian entanglement of formation 
of all three-mode mixed Gaussian states that are reductions of 
four-mode pure Gaussian states. Our results unveil the gen- 
eral structure and nature of correlations in continuous vari- 
able Gaussian states, and can find immediate applications in 
assessing the performance of quantum information primitives. 

Quantum discord |7| originates from the discrepancy be- 
tween two classically equivalent definitions of mutual infor- 
mation, a measure of total coiTelations in a quantum state. For 
classical probability distributions, the quantities I{A : B) = 
H{A) + H{B) - H{A, B), J{A : B) = H{A) - H{A\B), 
and J{B : A) = H{B) - H{B\A) all coincide due to 
Bayes' rule, where H is the Shannon entropy and the condi- 
tional entropy H{A\B) is an average of the Shannon entropies 
of A conditioned on the alternatives of B. For a bipartite 
quantum state qab, the mutual information can be defined as 
I{qab) = S{qa) + S{qb) — S{qab), where S stands for the 
Von Neumann entropy, S{g) — — Tr {glog g) (throughout the 
paper, log denotes the natural logarithm). The quantum ana- 
logue of J{A : B), known as one-way classical correlation 
and denoted as J^{gAB), is operationally associated with the 
distillable common randomness between the two parties 1 13 1, 
and depends on the measurements {11;}, 11^ = 1, made 
on i? 1 10|. The state of A after the measurement is given by 

Qa\i = B{QAB^i)/Pi, Pi = T^'c A,B{9AB^i)- A quantum 
analogue of the conditional entropy can then be defined as 
'H\Yii}{A\B) = '}2nPiS{gA\i): and the one-way classical cor- 
relation, maximized over all possible measurements, takes the 
form J^{gAB) = S'(£<a) - inf{ni} "Hjn J (^l-B). The quan- 
tum discord is finally defined as total minus classical correla- 
tions: 

'D^{QAB)^I{gAB)-J^{gAB) (1) 

= S{gB) - S{gAB) + mf H{nj(^|S)- 

We denote by J^{gAB) and T>^{gAB) the corresponding 
(generally different) quantities where the roles of A and B 
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are swapped. On pure states, quantum discord coincides with 
the entropy of entanglement S{gB) = S{gA)- States with 
zero discord represent essentially a classical probability dis- 
tribution embedded in a quantum system, while a positive 
discord, even on separable (mixed) states, is an indicator of 
quantumness, that arises e.g. when qab has entangled eigen- 
vectors 161 [m, and may operationally be associated to the 
impossibility of local broadcasting |4 |. 

We now define a Gaussian version of quantum discord 
and calculate it analytically for all two-mode Gaussian 
states. A two-mode Gaussian state qab is fully specified, 
up to local displacements, by its covariance matrix (CM) 
(TAB of elements aij — Tr[gAB{Ri, Rj}+] where R = 
{xa,PAiS:b,Pb) is the vector of phase-space operators sat- 
isfying the canonical commutation relations [Ri, Rj] = iflij, 

with n being the symplectic matrix (_°|^ q)*^ EJ. By means 
of local unitary (symplectic at the CM level) operations, ev- 
ery two-mode CM can be transformed in a standard form with 
diagonal subblocks 
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with 



a = al,l3 = bl 
7 = diag{c, d} 



(2) 



Let us define the symplectic invariants A = det a, B ^ 
det (3, C — det 7, and D = det ctab- The CM corresponds 
to a physical state iff A^B > \ and v± > 1, where the sym- 
plectic eigenvalues are defined by 2v^ = A ± ^/ /S? — AD 
with A = A + -B + 2C. A Gaussian state with CM ctab is 
entangled iff < 1, where the smallest symplectic eigen- 
value i>- of the partially transposed CM is obtained from 
by replacing C with — C, i.e. by time reversal 1 15 1. 

Both the one-way classical correlations J'^ and the quan- 
tum discord are entropic quantities, and therefore invari- 
ant under local unitaries. Hence, we can derive their closed 
formulae exploiting the standard form of a general two-mode 
Gaussian state, and later recast our results in terms of the 
four invariants of the state. The Gaussian quantum discord 
of a two-mode Gaussian state qab can be defined as the 
quantum discord where the conditional entropy is restricted 
to generalized Gaussain POVMs on B. These are all mea- 
surements executable using linear optics and homodyne de- 
tection He). We then have 'D^[qab) = S{qb) ~ S{gAB) + 
miiig(^) J d'ripB{ri)S(QA^)- Here the Gaussian measurement 
n_B on subsystem B can be written in general as Hb {rj) = 
Tr~^WB{r])'n°B^Biv) where WBiv) = e:>q){r]b^ -ri*b) is the 
Weyl operator, b = (xb + ipB)/V^, 7r~^ / d'^riIlB{ri) = 1 
and is the density matrix of a (generally mixed) single- 
mode Gaussian state. The conditional entropy is a concave 
function of the POVM elements, i.e. it is concave on the set of 
single-mode Gaussian states 11^. Gaussian states do not form 
a convex set, yet every Gaussian state admits a convex decom- 
position into pure Gaussian states, so it sufficient (as in the fi- 
nite dimensional case) to restrict to states 11^ that are pure, 
single-mode Gaussian states ifTTl whose CM we denote as 
ctq. The conditional state of subsystem A after the mea- 
surement IIb (r/) on B has a CM independent of the measure- 



ment outcome fT6l and given by e = a — j{P + ctq)" 
Recalling then that the Von Neumann entropy of a n-mode 
Gaussian state with CM cr can be computed as ifTSl S{(t) ~ 
J2i=i fi^i)' where i/i are the symplectic eigenvalues of the 
state and fix) = (^) log [3^] - {^) log [^], the 
one-way classical correlation and the Gaussian quantum dis- 
cord for two-mode Gaussian states with CM (tab 

•D^iTAB) - fiVB) - f{u^) - f{u+) + inf /(Tdete) . 

J^{(tab)^ f{^/A)~mif{VA^e). (3) 

To get closed formulae we need to minimize dct(£) over all 
CMs (To corresponding to pure one-mode Gaussian states, 
i.e. rotated squeezed states: ctq — i?(0)diag{A, 1/ X]Bj'^ [9), 
where A > and R(e) = (_™f,feto!e) ■ ^ general two- 
mode Gaussian state in standard form, one has E{X,9) = 
det(e) = [2a^{b+X){l+bX)-a (c^ + d^) {2b\ + + l) + 
a (c2 - d^) (A^ - 1) cos(20)+2c2d2A]/[2(6 + A)(l + bX)]. 
We set without loss of generality c > \d\. We now look for 
stationary points of E by studying its partial derivatives. One 
such point is A = 1, for which ctq is the identity (regardless 
of 6), and is a saddle point except when d — ±c. Next, the 
equation d\E — is quadratic in A, but one of its roots is 
always negative. The other root, given by A = A2 = [a6((i^ — 
c2) + c\d\ y/{a - a62 _^ 5^2) (a - ab^ + bd^)] / [ab^c^ - (a + 
bc^)d'^] and 6* = (or, equivalently, I/A2 and 6 = 7r/2), 
is a local minimum and is acceptable provided that A2 > 0, 
that is when 5 = {-acl +b [ab — cf) c|) > 0. Additional 
candidates for inf E have to be sought at the boundaries of 
the parameter space: a potential minimum lies at A — t- 0, 
6 ^ (01 equivalently A 00, 9 = 7r/2). In the whole 
physically allowed region for the CM parameters a, b, c and 
d, we have £'(1,0) > E{X2,0), with equality holding only 
when d ~ ±c, and i?(A2, 0) < E{0, 0). Thus, for any (Tab, 
info.„ det(£) is equal to £'(A2, 0) if 5 > and to E{0, 0) oth- 
erwise. In terms of symplectic invariants, the Gaussian quan- 
tum discord and the one-way classical correlation for a general 
two-mode Gaussian state ctab are given by Eq. Q with 

(4) 



E" 



inf det(e) 



2C^ + {-i+B){-A+D)+2\C\y/C^ + i-l+B)(-A+D) 



i-l+By 

AB) 2 < (1 + B) C2 (A 



D) 



AB-C^+D-y/C*+{-AB+Dy-2C^(AB+D) 

2B ' 

otherwise. 

Notice that this only depends on |C|, i.e., entangled (C < 0) 
and separable states are treated on equal footing. For states 
falling in the second case of Eq. Q, homodyne measurements 
(projections onto infinitely squeezed states, A = 0) on i? min- 
imize the conditional entropy of A. An example is when 



A = D 



C = (1 -B)/2, B^b^, 



(5) 



with 1 < 6 < 2a — 1, which is a mixed state of partial mini- 
mum uncertainty, i.e., one of its normal modes is the vacuum: 
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//_ = 1 f9l. On the other hand, the first case corresponds 
to a more general measurement, i.e., projection of mode B 
onto a squeezed state with unbalanced, finite variances on xb 
and pb- A notable class of states satisfying the first case are 
squeezed thermal states (including pure states), characterized 
by d = ±c, for which the conditional entropy is in particu- 
lar minimized by heterodyne measurements (projection onto 
coherent states, A ~ 1). In general, Gaussian quantum dis- 
cord can be accessed experimentally by linear optics, and our 
finding provides the optimal measurements to verify quantum 
correlations given the CM of a two-mode Gaussian state. 

We now analyze the relationships between classical corre- 
lations, quantum discord, separability and entanglement. For 
every entangled state the quantum discord is strictly posi- 
tive (since S{gB) ~ S{gAB) > 0). Almost all separable 
states in finite dimensions have also nonzero discord |5|. In 
any dimension (including infinite dimensions under the con- 
straint of finite mean energy), the states qab with zero dis- 
cord are the ones that saturate the strong subadditivity inequal- 
ity for the Von Neumann entropy on a tripartite state qabc 
where C is an ancillary system realizing the measurements on 
B llT4l |T9| . From the characterization of such states in the 
Gaussian scenario [ 20] (see Supplementary Appendix A ETI 
for more details), it follows that the only two-mode Gaussian 
states with zero Gaussian quantum discord are product states 
(Tab = q; ® /3, i.e., states with no correlations at all, that 
constitute a zero measure set. Quite remarkably, then, all cor- 
related two-mode Gaussian states have non-classical corre- 
lations certified by a nonzero quantum discord. This is in 
qualitative agreement with a recent study | |22| where it has 
been demonstrated numerically, using the lack of a positive- 
definite Glauber-Sudarshan P representation as a nonclassi- 
cality criterion, that essentially all two-mode Gaussian states 
are a priori nonclassical. 

For Gaussian states with asymptotically diverging mean en- 
ergy, however, interesting correlation structures arise. Con- 
sider the squeezed thermal state given by 

a = cosh(2s), h = cosh^ rcosh(2s) + sinh^ r, 
c = —d = coshr sinh(2s) . 

For r = 0, this is a pure two-mode squeezed vac- 
uum state, whose entanglement is an increasing function of 
s. In the limit r — > oo, it is asymptotically separable 
(but not in product form). Concerning the discord (mini- 
mized in this example by heterodyne detections), we find 
'D^{(Tab) = /[cosh^ r cosh(2s) + sinh^ r] — /[cosher + 
cosh(2s) sinhV] ' > and X>^(cr^B) = /[cosh(2s)]- 
/[cosh^ r+cosh(2s) sinh^ r]+/[cosh(2r)] ''''^°°> 1. While 
these limiting values are associated with ideal, unnormaliz- 
able states, they can be approached arbitrarily close by physi- 
cal Gaussian states with large, but finite mean energy. Hence, 
surprisingly, there exist bipartite Gaussian states such that: (i) 
they are non-product states, with arbitrarily large correlation 
matrix 7, yet have infinitesimal quantum discord; (ii) their 
quantum correlations can be revealed by probing only one 



subsystem, but not the other Thus motivated, we have ex- 
plored the discord asymmetry for one million randomly gen- 
erated (separable and entangled) two-mode Gaussian states. 
Let V''''' = maxjl?^,!?^}, X>'"'" = min{I?^,2?^} for 
a given CM. We find numerically that P'"^'' - P™''^ < 
X>™'"/[exp(X'""") - 1] < 1. The leftmost bound is saturated 
by states of Eq. (j6]l in the limit s — > 00, and unity is reached 
for r — >■ 00 as well. The maximum discord asymmetry decays 
exponentially with I?™'", so when the discord calculated in 
either way is large, we have de facto T)^ = T)^ . 

Next we ask: To what extent can separable Gaussian states 
be quantumly correlated? While their discord is typically 
nonzero (but for product states), we find that it cannot ex- 
ceed one unit of information [Fig.[T|left)]. In the Supplemen- 
tary Appendix B II2TI we prove that for all two-mode sep- 
arable Gaussian states, I?^((t^^]^) < [{b — l)/2]log[(fe + 
l)/(^ ~ 1)] — 1- The first inequality is saturated by separable 
squeezed thermal states whose correlation matrix has maxi- 
mum determinant C and whose CM has maximum asymmetry 
between the two modes: c — d~ 1+ab—a—b, a — 00 [solid 
(red) curve in Fig. [TJleft)]. The second bound is reached for 
& — 00. This result immediately implies a sufficient condition 
for the entanglement of Gaussian states given their discord: ;/ 
V^^ctab) > 1, then (tab is entangled. 

We now focus on entangled states, and study how Gaussian 
quantum discord compares quantitatively to the entanglement 
of the state, specifically measured by the most "compatible" 
measure available, the Gaussian entanglement of formation 
(Gaussian EoF) £g llli This is defined for Gaussian states 
QAB as the convex roof of the Von Neumann entropy of entan- 
glement, restricted to decompositions of qab into pure Gaus- 
sian states. It can be evaluated via a minimization over CMs: 

£g{o-ab) = info-^^<crAB : dot(,T^^)=i/(\/deta')' where 
the infimum runs over all pure bipartite Gaussian states with 
CM ct'ab smaller than (tab, and a' is the reduction of (t'ab 
corresponding to the marginal state of mode A. Compact for- 
mulae for £g exist for all symmetric two-mode states (where 
the Gaussian EoF coincides with the true EoF as the Gaus- 
sian decomposition is optimal) 123)1 . as well as asymmetric 
ones ill] El. hi Fig. [TJright), we plot vs for 30000 
randomly generated two-mode Gaussian states. We find that 
for a given entanglement degree, the discord is bounded both 
from above and below. To find the upper bound analytically, 
we can restrict, as in the separable case, to squeezed thermal 
states (see ET]) with d = — c. Further optimization within 
this family of states yields that, for all two-mode entangled 
Gaussian states, the quantum discord satisfies: 

V^icTAB) < max{£G('7'AB), 2cosh^ r log(cothr)} , (7) 

where £g(cab) = /(I + 2sinh~^r) implicitly defines 
r. The rightmost bound [solid (red) curve in Fig. [Tfright)] 
dominates in the low entanglement regime (£c < 2 log 2) 
and corresponds to of the states of Eq. (|6]l in the limit 
s ^ 00. The leftmost bound in Eq. (|7]i [dotted (green) 
line in Fig. [TJright)] is instead reached on pure states, and 
sets an upper limit to the quantum discord of all two-mode 
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separable states entangled states 




FIG. 1. (Color online). Left: Gaussian quantum discord versus 
marginal entropy for separable two-mode Gaussian states. Right: 
Gaussian quantum discord versus Gaussian EoF for entangled two- 
mode Gaussian states. See text for details of the bounding curves. 

Gaussian states with sufficiently high entanglement. On the 
other hand, for a given £g, the Gaussian discord also satis- 
fies 2?^((t^b) > 21og(cothr), with r as before. This lower 
bound follows from the fact that states of Eq. (|6]l with s ^ oo 
are extremal for the discord asymmetry, and corresponds to 
of those states [dashed (blue) line in Fig. [TJright)]. In- 
terestingly this entails that, asymptotically, for all two-mode 
Gaussian states with £q 3> 0, their Gaussian quantum dis- 
cord lies between Eq — 1 and Eq. 

A further key result of our study is that Eq. (j4|i provides a 
closed, computable formula for the Gaussian EoF of a class 
of three-mode mixed Gaussian states. This is possible thanks 
to a duality relation between (Gaussian) classical coiTela- 
tions and (Gaussian) EoF lfT2l (see also ll25l ). Let qabst 
be a purification of qab, i-e. a pure (Gaussian) state such 
that Ti:st\qabst\ = Qab (we need in general two ancil- 
lary modes S and T to construct such a purification ifTSl ). 
Then: J^{qab) + ^{qast) = S{qa) where E{qast) de- 
notes the EoF between party A and the block of modes ST. 
In the Gaussian framework, from Eq. (j3]l we have simply: 
Eq{(tast) — info-o /(Vdet e). The states with CM ctast 
encompass all three-mode Gaussian states that are reductions 
of a four-mode pure Gaussian state. Their symplectic spec- 
trum is of the form {1,1,6}, i.e. they are mixed states of par- 
tial minimum uncertainty, with two vacua as normal modes. 
For all such states, we now present an analytic method to 
compute the Gaussian EoF across the bipaitition A x {ST): 
first, construct a purification, i.e. append an ancillary mode B 
(with det(<TB) = det(<T^sT) = fc^) such that (Tabst is pure 
ifrniT? !. Then, evaluate £;™'" of the marginal state ctab from 
Eq. (Wji. Finally, /(V £'™'") is the Gaussian EoF between A 
and ST. An example of such a state, of relevance in a crypto- 
graphic setting, is provided in fl\\ (Appendix C). 

Finally, we wish to point out that for states ctab with 
~ 1, the purification requires a single ancillary mode S, 
and the Gaussian EoF between modes A and S, as computed 
through Eq. agrees with the formula derived in |24|. In 
the special case of these states given by Eq. (jSjl, the comple- 
mentary state of modes AS is symmetric. Thus the optimal- 
ity of Gaussian decompositions for the EoF L23 1 implies that 
no non-Gaussian measurement can do better than the optimal 
Gaussian POVM devised here for the optimization of the con- 
ditional entropy. In other words the (unrestricted) quantum 



discord and one-way classical correlations in the family of 
two-mode states (Tab of Eq. (Q coincide with their Gaus- 
sian counterparts. For general two-mode Gaussian states, 
it is an open question whether non-Gaussian measurements 
(e.g. photo-detection) can lead to a fuither minimization of 
the discord. Insights on this issue might be drawn if the addi- 
tivity conjecture forbosonic channels were proven |25|. 

This Letter opens the way for the study of quantum discord 
in general multimode correlated bosonic systems, and through 
the paradigmatic two-mode case demonstrates the 'truly quan- 
tum' nature of Gaussian states, reinforcing their key role in 
continuous variable quantum information processing. Such 
studies could also lead to a deeper understanding of entangle- 
ment and the general structure of distributed correlations in 
harmonic lattices. 
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Note added. — After the completion of this study, during the 
writing of the present paper, another work appeared ll26l 
where Gaussian quantum discord is independently defined for 
Gaussian states and explicitly calculated only in the specific 
case of two-mode squeezed thermal states. 
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A. Gaussian states saturating the strong subadditivity of entropy 

For a general state qabc of a tripartite system, the Von 
Neumann entropy is strongly subadditive: 

S{qab) + S{qbc) > SiQABc) + S{gB) ■ 

For finite-dimensional and infinite-dimensional systems (in 
the latter case, under the restriction of finite mean energy), it 
can be proven that the states qab with zero quantum discord 
V are obtained from tripartite states qabc Qab ® Qc that 
saturate the strong subadditivity. Here C is an ancillary sys- 
tem needed to realize the most general measurement on sub- 
system B, and c[qabc\ = Qab- lfT4l[T9l . In the Gaussian 
scenario, C is in general a collection of Gaussian modes, and 
interactions are restricted to Gaussian ones, i.e., symplectic 
operations and homodyne detection, that preserve the Gaus- 
sian character of the involved states. Therefore, the Gaussian 
states with CM aAB that have zero Gaussian quantum dis- 
cord are the ones such that strong subadditivity is saturated on 
Gaussian extensions otabc 7^ (^ab ffl cc- In the general case 
of tripartite Gaussian states where each subsystem A, B, and 
C contains an arbitrary number of modes, Petz and Pitrik have 
recently characterized the subset of Gaussian states (tagged as 
Markov states) saturating the strong subadditivity EOll . For a 
tripartite n-mode CM 

fABC = (^ ai, . . . , ai , ai+i, . . . am , am+i , ■■■an 

A B C 

corresponds to a Gaussian Markov state if and only if there 
exists a splitting within the block of modes B (including the 
boundaries) such that the CM is block-diagonal, i.e., if and 
only if 

f ABC = f^ai,...,ai^j ® f ai+j^i ,. . . ,a„ 

for some j, < j < m — I. In our setting, the block B 
is made of a single mode only, hence the only possibilities 
are j = 0,1. Consequently, the only Gaussian Markov states 
in this case are either of the form ctabc — cta ffi ctbc or 
CABC = CAB ffi cTc- The latter having been discounted al- 
ready, it follows that the only zero-discord two-mode Gaus- 
sian states (under the assumption of finite mean energy) are 
the reductions of (ta(B(Tbc after partial trace on C, i.e., prod- 
uct states (Tab = ca © ""b- 



B. Upper bound on the quantum discord of separable Gaussian 
states at given marginal entropy 

Here we prove the upper bound on quantum discord for sep- 
arable Gaussian states at fixed B = b^. First, we observe that, 
since the discord 'D^{crAB) = VIA, B,C, D) (notice the 
slimmed notation and the explicit dependence on the invari- 
ants) involves a minimization over the Gaussian POVMs on 
system B, any such measurement provides an upper bound 
for the discord. In particular, heterodyne measurements do 
the job, being V{A, B, C, D) < Vhetero{A, B, C, D) for all 
two-mode Gaussian states. Since we know that squeezed 
thermal states are the only ones that saturate the bound, we 
can restrict to them for the maximization of the discord, and 
write that for fixed A, B, and C, it is V{A, B, C, D) < 
ViA, B, C, {y/AB± C)2). The ± sign distinguishes between 
the two cases, d = =Fc, that can be associated to separa- 
ble squeezed thermal states. Since the conditional entropy, 
from Eq. Q, does not depend on the sign of C, the maxi- 
mal discord is obtained in the case of minimal global entropy 
S{crAB)^ From the expression of the symplectic eigenval- 
ues iy± and the form of the concave function f{x) we find 
that the case d = c, i.e. D = {\/ AB — C)^ maximizes the 
global entropy. Next, we observe that the discord for the fam- 
ily of separable squeezed thermal states under investigation is 
a monotonically increasing function of C, hence it is maxi- 
mized at the boundary C = 1 — ^/A — \/B + \J AB. We have 
thus proven that, for all two-mode separable Gaussian states, 
V{A, B, C, D) < V{a^,b^, 1 - a - b + ah, {a + b - l)^) = 
f{b) -f{a + b-l) + /[(2a + b-l)/{l + b)]. This is by itself 
an useful upper bound that depends on both local entropies of 
the two-mode state. To get a function of b only, we observe 
that the obtained bound is a monotonically increasing function 
of a, and is hence maximized at a oo. Using the fact that 
f{y) ~ f{^) logijj/x), we finally find the bound in 

the text: V^{(tYb) < [ib-l)/2] log[(fe + l)/(6 - 1)] which 
is an increasing function of 6 converging to 1 for 6 ^ oo. 

C. Example of a three-mode mixed Gaussian state whose 
Gaussian EoF can be computed exactly 

We now provide an example of a three-mode family of 
states with CM <tast, of relevance in a cryptographic set- 
ting, whose normal mode decomposition accommodates two 
vacua. Suppose that Sally and Tom are sharing a pure two- 
mode squeezed state of modes S and T, cr'glp, of the form 
Eq. ^ with a ^ b = s,c — —d = \/ — 1, as a quan- 
tum communication channel. A malicious Eve remotely pre- 
pares a two-mode squeezed state c™^ of her modes A and 
B, given by a = 6 = ri, c = —d = ^/-n? — 1. Eve then, 
maintaining local access to her mode B, attempts to inter- 
cept the communication by sending her mode A to interfere 
with Sally's mode S through a beam-splitter with transmit- 
tivity t, as schematically depicted in Fig. |2] After the at- 
tack, the global state of the four involved modes has a CM 
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FIG. 2. Schematic picture of the production of the four-mode 
(tabst discussed in the text. Sally and Tom control modes S and 
T, respectively, while Eve owns modes A and B. 



given by (tabst = BAs{t){(T'J^B ® (Tfrj.)BAsit), where 
Bij {t) can be expressed in the standard form of Eq. (j2ji with 

a = —b = ^/i, c = d = \/\ — t. Now Eve's ability to 



eavesdrop depends crucially on the entanglement between her 
probe mode A and the target modes ST. The reduced state 
(TAST has symplectic spectrum {1, 1, n} so it belongs to the 
class whose Gaussian EoF can be calculated thanks to our re- 
sult. We find that £g{<^ast) ^ f[t + o,{1 — t)]. This is equal 
to the conditional entropy of the probe mode A after Eve per- 
forms a heterodyne measurement (optimal in this case) on her 
local mode B. In other words, Eve generates quantum correla- 
tions with the communicating parties, no matter how weak the 
entanglement (quantified by n) she had established between 
her modes AB in the first place. For t ^ 0, there is an entan- 
glement swapping between Eve's modes and Sally and Tom's 
modes, yielding the most destructive attack. In the limit of 
two perfect quantum channels, a,n oo. Eve is sacrificing 
the (infinite) entanglement between her two modes, spending 
it to get (infinite) entanglement with the two spied parties: yet, 
some quantum correlations persist between A and B, being 
'D^{crAB) = 1 — log 2 in this limit. 



